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Geometric phases and information flows of a two-level system coupled to its environment are 
calculated and analyzed. The information flow is defined as a cumulant of changes in trace distance 
between two quantum states, which is similar to the measure for non-Markovianity given by Breuer. 
We obtain an analytic relation between the geometric phase and the information flow for pure 
initial states, and a numerical result for mixed initial states. The geometric phase behaves differently 
depending on whether there are information flows back to the two-level system from its environment. 
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I. INTRODUCTION 



Geometric phase has not been noticed for a long time until Pancharatnam's study [1} and Berry's discovery 0. 
Shortly, Simon gave a geometric interpretation of this kind of phase in the language of differential geometry and fibre 
bundles Q • Since then there was a keen interest in holonomy effect in quantum theory, which leads to many extensions 
of the geometric phase, including the geometric phase acquired in a non-adiabatic and cyclic evolution by Aharonov 
and Anandan Q, in a non-adiabatic and non-cyclic evolution by Samel and Bhandari Q and in a non-adiabatic, 
non-cyclic and non- unitary evolution by Mukunda |6]. All those investigations were focused on pure state. For a 
practical quantum system, however, its state would be mixed due to the unavoidable coupling to its environment. 
^/ This motivates the study on the geometric phase for mixed states [7j-[Toj ■ which was defined by [lOj . 
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where £j(t) and \ipi(t)) are the eigenvalues and the corresponding eigenstates of the density matrix, respectively. 
Throughout this paper, we will use this definition to study the geometric phase. 

On the other hand, it is complicated to exactly describe the dynamics of open systems due to the huge number 
of variables in an environment. In order to simplify the description, one could make some approximations such as 
the weakly coupling and the Markovian approximation. Under these approximations, we can obtain a Markovian 
master equation that describes the dynamics of the open system without memory of its history. However, many 
systems exhibit strong non-Markovian effect and can not be described by the Markovian master equation. It is 
then interesting to study the geometric phase in a non-Markovian dynamics, and establish the relation between the 
geometric phase and the non-Markovianity of the dynamics. 

The non-Markovianity may be defined in many ways [34I3] j for instance, in Ref. [l4[ the authors proposed a scheme 
to quantify the degree of the non-Markovianity based on the trace distance of two quantum states (T^J, and in Ref.fTEj 
the non-Markovianity was measured by exploiting the specific traits of quantum correlations. The effects of non- 
Markovianity on geometric phase has been considered by several works, e.g. Refs. [12j and [13j . In this paper, we will 
establish a relation between the geometric phase and the information flows. We divide the information flow into two 
types, i.e., the information flow from the open system into the environment M. (call forward information flow) and the 
information flow back from the environment to the open system M (call backward information flow) . The backward 
information flow M in fact is a modified measure of non-Markovianity given by Breuer and his co-workers. This paper 
is organized as follows. In Sec. II, we present a definition for the information flow based on the trace distance, then we 
establish the relation between the information flow and the geometric phase for pure initial states. The case of mixed 
initial states is considered in Sec. Ill, where the geometric phase and the information flow are calculated numerically. 
Finally, we present conclusion and discussions in Sec. IV. 
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II. INFORMATION FLOW AND GEOMETRIC PHASE FOR PURE INITIAL STATES 

A. A measure for information flow 

Here we first recall the measure for non-Mar kovianity denned by Breuer[3|. This definition is based on the so-called 
trace distance between two states p\ and p 2 

D(pi,P2) = lTr\\ Pl - P2 \\, (2) 

where \ \A\\ — \J AA^ . For a two-level system, this trace distance is equal to one half of the ordinary Euclidean distance 
between the two states on the Bloch sphere, i.e. D(pi, p 2 ) = h\fi — r 2 \, where fj is the Bloch vector for state pj. The 
change rate of the trace distance can be represented as 

a( Pl (t), P2 (t)) = j t D( Pl (t),p 2 (t)). (3) 

When a < 0, pi(t) and p 2 (t) approach to each other in the dynamics and this can be understood as information flow 
from the system to the environment; when a > 0, pi(t) and p 2 (t) is away from each other in the time evolution, 
and this can be interpreted as information flow back to the system, which is treated as a typical character of the 
non-Markovianity. As shown in Refs.[l4[ and [l8j |. one can define a measure of non-Markovianity Mb (t) for a quantum 
process *f?(t) by maximizing over all initial states (pi(0), p 2 (0)) of total gain of the trace distance, namely, 

A/" B (1<) = max / a{ Pl {t),p 2 (t))dt. (4) 

Pi,a(0) J a>Q 

The measure of non-Markovianity only characterizes the gain of the trace distance in the dynamics. To describe the 
loss and gain of the trance distance uniformly, we extend the concept of information flow in the following. In fact, one 
of the information flows defined below is a modified measure of non-Markovianity in Eq.(j4j. We choose a steady state 
as p 2 in Eq. ([2]) and call it standard state. This choice makes the maximization easy, however, it can not measure all 
non-Markovian dynamics. Fortunately, this simplified measure of non-Markovianity is in agreement with the measure 
given in [14j in our case. 

Now we define M. and M to measure the information gain and lose for a given initial state p\ (0) in the dynamics 
^(t). We will refer this information gain and loss uniformly as formation flows. 

M(V) = I a( Pl (t),p 2 )dt, (5) 

JiT>0 

and 

M(*) = -[ a( Pl (t),p 2 )dt. (6) 

Obviously, the above two measures satisfy, 

D{ Pl {t),p 2 ) = D( Pl (0),p 2 ) +M(p 1 (t),p 2 )-M(p 1 (t),p 2 ) . (7) 

The difference between M(t) and Mb is as follows. Mb is the maximum information flow back to the system in the 
dynamics. Hence it does not depend on the initial condition. However M{t) is the information flow back to system 
with respect to the standard state p 2 {Q) in the interval (0, t) (assumed that a steady state for quantum process 
exists, and the standard state p 2 (0) is exactly the steady state as before). This simplification is true in the situation 
considered in this paper, i.e., a two- level system coupling to vacuum electromagnetic fields at zero temperature with 
standard state p 2 in the Bloch sphere representation as r 2 = (0,0, —1). 

With these definition and notations, we will discuss the relationship between the geometric phase and the informa- 
tion flows between the open system and its environment. 

B. The geometric phase for pure initial states 

Consider a two-level system coupled to its environment at zero temperature. The general form of the density 
matrix can be expressed as p(t) = i(l + f(t) ■ <r), where o = (<t x ,<j v ,o z ) is the Pauli matrices, and r(t) = r(t) ■ 
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(sin#(t) cos</>(i), sin#(i) sin <j>(t), cos9(t)) is the Bloch vector. For pure initial states |r(0)| 2 = 1, while for mixed initial 
states |f(0) | 2 < 1. It is easy to obtain the instantaneous eigenvalues of the above density matrix p(t) as, 

e±(t) = ~(l±Kt)|). (8) 

It is obvious that for the pure initial states, the eigenvalue e_ (i = 0) = 0, which means that the eigenstate correspond- 
ing to the eigenvalue £-(i) gives no contribution to the geometric phase. This simplifies our calculation and helps us 
to obtain an analytic result for the geometric phase. The eigenstates corresponding to the eigenvalues in Eq.© can 
be written as, 

!*+<«»- {J$$m ) •"*-«»- ) • » 

Now substitute Eq.© and Eq.© into the Eq.([T]) with an assumption that <j) = cj £ + <j>o (where ojq and 0o are 
constants, this is reasonable for different kinds of master equation [l2|, EH), the geometric phase with pure initial 
states (their Bloch vector is r(0)) can be obtained as 

$gp = - / w cos 2 ^dt, (10) 
Jo z 

where we set T — (n = 1, 2, 3, ...), i.e., a multiple of the quasi-period 

Next, we establish the relationship between the geometric phase and the information flows M.(t) and M{t). Keeping 
the relationship Eq. ([7j) between the N{t) and M. (t) in mind, we obtain the geometric phase acquired by the system, 



$ GP =- / uo[-+ — rz(t) =]dt. (11) 

Jo 2 2^fA[D{Q)+N{t)-M{t)] 2 -2r z {t)-r 

where r z (t) — r{t) cos9(t) is the z-component of r(t) and D(0) = D(pi(0), pi). It is shown that when the system is 
closed, i.e. cosd(t) = cos#o and r(t) — 1, where #o is initial polar angle on the Bloch sphere, and setting T = 2w/uio, 
Ea. (fni) reduced to the well-known form, $^p 0sed - ) = — 7r(l + cos^o) 5 that is the geometric phase acquired by a two-level 
quantum system in a rotating magnetic field. For a Markovian process, Af(t) is always zero and Ai(t) increases with 
time until it approaches D(0). In this case, the geometric phase is only influenced by the information flow to the 
environment. When we consider the non-Markovian effects, the situation is more complicated, and the information 
flow Af(t) back to the open system has a link to the geometric phase given by Eq. ([TTj) . We will discuss it numerically 
in Sec.III. 



III. THE GEOMETRIC PHASE FOR MIXED STATE 



In this section, we will study the geometric phase of a two-level system with mixed initial states. Because it is 
difficult to get an analytical result like Eq.(|lip for the information flows and the geometric phase, we here numerically 
and perturbatively establish a relation between the geometric phase and the information flow. The perturbation is 
carried out to the first order in the coupling constant, while the numerical results are for a wide range of coupling 
constants. Two types of master equation, the time-local master equation and the memory kernel master equation 
with exponential memory, will be considered. 



A. The time-local master equation 



Let us consider a two-level system interacting with a vacuum field at zero temperature whose spectral density is 
Lorentzian [HE [11, 

1 W 2 X 

7r (wo — ^) + * 

Here W is the coupling constant between the system and the environment, ujq is the atomic transition frequency 
which is of the time scale tq ~ w^ 1 , and A is the spectral width of the coupling that is connected to the environment 
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correlation time, tb ~ A 1 . The dynamics of this system is governed by the following master equation (namely the 
time-local master equation) [Til ]. 

p{t) = -iA(t)[a + a^,p(t)} + T{t)[2a.p{t)a + - tr+a-p(t) - p{t)a + a.}, (13) 

where a± are the pauli operators and the parameters 

A(t) = -$[§§], r(t) = -3?[|g], (14) 

play the role of Lamb shift and decay rate for the system, respectively. Here c(t) can be calculated by means of 
Laplace transform as, 



c(t) = exp 



(A + icoo)t 



, fit a . ot\ 

cosh y + - smh — ) , (15) 



with f2 = VA 2 — 4VK 2 . We note that the rate R = ^ indicates the strength of the non-Markovianity. With R < A 
the dynamics is called time-dependent Markovian, while for R > i the dynamics is non-Markovian. Here we assume 
again that the initial state of the open system is 

p(0) = ±(l+f -<7), (16) 

where ro = ro • (sin#o cos </>o, sin#o sin </>o, cos^o) and |rb| 2 < 1. With this initial condition, the density matrix of the 
system at time t can be obtained from the master equation Eq. p^|) 

, , _ 1 / (1 + r cos6» )|c(t)| 2 r o sin0oexp(#o)c(t) \ 
P[> ~ 2 \rosm0 o exp(-i<h)c*(t) 2 - (1 + r cos 8 )\c(t)\ 2 J ' 

The eigenvalues and the eigenstates of the reduced density matrix Eq. (IT71) can be easily obtained as 

e±(i) = 5(l±r(t)), 
llM*)> = (^sSX**)). ( 18 ) 



-cosf. 



where tan&> t = and r(i) = ^[(i + rQ cos 9 Q )\c(t)\ 2 I] 2 + r 2 sin 2 o |c(t)| 2 . 

We can expand the geometric phase with respect to the coupling strength W 2 up to the first order (i.e., in the weak 
coupling limit), that is 

$ GP (T) = $gp - ^ 2 [tan$^C 1 (r o ,0 o ) Kl (A,T) + Wo cos- 2 $^C 2 (r o ,0 o ) K2 (A,T)], (19) 
where Kl (A,T) = d -^f \ w . =0 = i~°*p(-*H _ T = ^ aj^^ ^ = t + ^[ exp( _Ar) - 1] - 

rp2 (Q\ t _ _ _ ...... (0) 

ov, and •I'qp is the geometric phase acquired under the unitary evolution with mixed initial states, $q P = 
arctan(ro tan(— m{l + cosOq))). The parameter Ci(ro,8o)(i — 1,2) is a constant relative to the initial condition 
given by 

Ci(r ,6 ) = |(r + r cos 2 00 + 20036*0), 
C 2 (r ,do) = ±(l+ r ° sin2e 2 ° cose ° -cos 2 o ). (20) 

It is interesting to calculate the trace distance between p(t) defined in Eq. pT)) and the standard state in the weak 
coupling limit. Substituting Eq. (|17l) into Eq.Q, and expanding the trace distance up to the first order in W 2 , we 
have 

W 2 d I c(t) I 2 

D(t) = D(t)\ w *= + 4j, (f) i 2 _ q x [^(1 + cos 2 g ) + 27-0 008 00] ^ l^ 2 =°' ( 21 ) 
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FIG. 1: (Color online) The geometric phase (in units of 7r) as a function of the parameter R with different z. The coupling 
strength is W = O.luo and ujo = 1. The azimuthal angle of the two components of the initial state are chosen as $o = 7! "/4, 



where D(t) \w 2 =o= \Ao + 1 + ^ r o cos#o/2 and d g^)-j \w 2 =o= 1 cx p( At ) _ jL. D(i)| w/2=0 is the trace distance at 
time i under unitary evolution, it is not difficult to prove that D(t) \w 2 =o= D(0). Then according to Eq.([7|), we have 

W 2 3 I r(t) I 2 

M(t)-M(t)±—-—x [rg(l + cos 2 g ) + 2r Q cosg ] \ w2=Q . (22) 

If coupling strength is very weak, there is no information flow back into the system in a quasi-period, i.e. J\f(T) = 0. 
So it is straightforward to obtain, 

W 2 8 I r(t) I 2 

M(t) = - w ^-— x [r 2 (I + cos 2 do) + 2r cosg ] LUJ | W2=Q . (23 ) 

This result tells us that Ai(t) increases monotonically with the increasing of 1/A. Comparing this result with Ea. dT9t . 
one may find that the dependence of Ai(T) and &gp(T) on the spectral width A is almost the same for mixed initial 
states in the weak coupling limit. 

Numerical results for the geometric phase and the information flow M under the weak coupling limit is shown in 
Figffland FigfJJ The initial states are chosen as , 

p(0) = i^I + z|0(CI, (24) 

where |£) = cosi?o|0) + sin^o exp(zy>o)|l) is a pure state, z € [0, 1] and I is a 4 x 4 unitary matrix. For z — 0, the 
density matrix is the maximumlly mixed state, while they reduce to a pure one in the case of z = 1. In the language 
of Bloch vector, the initial state Eq. (|24|) can be represented as 

r = z,9 = 2tfo,<?f>o = ¥>o- (25) 

Now we establish the relation between the geometric phase and the information flow with different coupling strengths 
one by one. E.g., W = O.lwo, W = ojq and W = 10wo will be chosen to explore the geometric phase and the information 
flow, numerical results are shown in Fig. ([T])— (J5J) . Fig[T]and Fig[2]show the geometric phase $qp and the information 
flow M. as a function of R with different parameter z for weak coupling (W — O.Iwo). The geometric phase <&gp is 
plotted in units of w. We can see from Fig[T]that the geometric phase increases monotonically with the parameter R. 
In this case, the coupling strength W is small enough so that there is no information flowing back into the system in 
a quasi-period. For pure initial states, i.e. z = 1 , comparing FigQ]with Figf2l we find that the larger the geometric 
phase is, the more information flow to the environment, this is confirmed by Eq. ljlip . Because N{R) is always zero 
in a quasi-period, the geometric phase mainly depends on A4(R). 

With increase of the coupling strength W , the correlation time of the environment tb approaches to the time T, 
which indicates that Af(R) > in a quasi-period. In Fig|3] and FigJH we plot the information flow Af(R) and the 
geometric phase <&gp as a function of R for W — luq. Here we only focus on the information flow AT, because it 
characterizes the non-Markovianity of the open system and describe the backward information flow. 
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FIG. 2: (Color online) The information flow M as a function of the parameter R with different z. The parameters are the 
same as in FigfTJ 
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FIG. 3: (Color online) The information flow M as a function of the parameter R with different z. Here the coupling strength 
W = ljo and ujq = 1. The initial states of the open system are chosen as Fig[TJ 

Comparing Figj3]with Fig|4] we may find that, when the information flows back to the system, the geometric phase 
decreases with the increase of Af(R) and, in the region of Af(R) — 0, the behavior of the geometric phase is similar 
with the case in the weak coupling limit. This indicates the backward information flow (i.e., information flow back 
to the open system) affects the geometric phase acquired by the open system. This phenomena can be understood as 
follows: For this time-local master equation, when the information flows back to the system, the Bloch vector moves 
toward to the north pole of the sphere, then the geometrical phase which is interpreted as the solid angle in Bloch 
sphere decreases. 

When W is very large, the system will drop into steady state after time T. In this case, by comparing the geometric 
phase $gp an d the information flow J\f, we can see that, when Af(R) = (as shown in Figj6]), the behavior of the 
geometric phase (as shown in Fig[5]) is very similar to the case of W = 0.1cj and W = cj with Af(R) = 0. With 
the increase of R, the information flow back into the open system and the geometric phase decreases with i?, this is 
very similar to the case when W = u>o in the region of M(R) > 0. Based on these observations, we conclude that for 
a dynamics described by the time local master equation, if the geometrical phase is inversely proportional to R, the 
dynamics must be non-Markovian. In other words, the non-Markovianity can be reflected in the geometrical phase to 
a certain extent. This conclusion are both valid for pure and mixed initial states. 

It is seemingly that the point (in the R axis) where the geometric phase arrives its extremum is exactly the point 
where Af(T) begin to increase (in the following we will call it as the critical point), however, by carefully examination, 
we find that this is not the case. According to Eq. (fTTj) and the definition of trace distance, we can clarify that the 
critical point is the very point where the integrand of the geometric phase reaches its minimum (a detail of proof can 
be found in the Appendix). Moreover, the integrand of the geometrical phase behaves similarly with the information 
flowJV(t). 

The information flow JV(T) and the geometric phase $gp as a function of $o and R are plotted in FigEl When 
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FIG. 4: (Color online) The geometric phase as a function of the parameter R with different z. Here the coupling strength 
W — OJo and Wo = 1. The initial state of the open system are chosen as FigfT] 

J\f(T) > 0, it can be seen from the figure that the more the information flows back, the smaller the geometric phase; 
in the region of Af(T) = 0, no matter how to choose the initial azimuthal angle $o> the geometric phase increase 
with R, which is exactly the finding of our analytical analysis. When $o = 7r /2, both Af(T) and <I>gp reach its 
extremum, in this case, the geometric phase do not change with R and ^gp — 27r, but the information flow changes 
indeed. This can be understood as follow: when i?q = 7r /2, sin2^o = 0, the initial density matrix is diagonal, i.e. 
p(0) = diag{(l — z)/2, (1+z) /2}, with this diagonal density matrix, the geometric phase is always equal to 27r and M(T) 
changes with R since the information exchange between the system and the environment varies with R. The situation 
remains unchanged for the cases where #o = and d — tt. Furthermore we find that, although $gp = 27r for both 
i?o = tt/2 and do — tt, the information flows JV(T) are completely different. This can be explained as the difference 
in the initial states, which are p(0) = diag{(l - z)/2, (1 + z)/2} for i} = n/2 and p(0) = diag{(l + z)/2, (1 - z)/2} 
for -do = 7T. By the definition of the geometric phase, it depends on the spectrum of the density matrix, which are the 
same for the initial states, leading to the same geometric phase acquired in the dynamics. But the information flow 
Af(T) for $o = tt is larger than that for i?q — 7r /2, this is because the information flow was defined as the distance 
between the actual state and the standard state of the open system, which are different for the initial states. 




B. The memory kernel master equation with exponential memory 



Now we consider the geometric phase of an open two-level system governed by the memory kernel master equation 
with exponential memory. Here, we just apply this model to calculate the geometric phase but do not discuss the 
positivity of the master equation in detail. 

In the interaction picture, the integro-differential master equation with memory kernel can be expressed as 

p(t)= [ K(t')£p(t-t')dt' (26) 
Jo 

where C is Liouvillian superoperator which takes the form 

£-P = ^7o(2cr-pcr + - <7 + cr_p - po+<J-), (27) 

and 70 is the dissipation rate, K(t) represents the memory effect called Shabani-Lidar memory kernel (22l. [23| 

K(t) = 7exp(-7i). (28) 

We call tr = — the memory time. It is not difficult to solve this intcgro-diffcrcntial equation by the Laplace transform 
with the initial condition Eg. (1161) . In the Schrodinger picture, the solution is 



Pit) ^ 



(1 + ro cos0 o K(C,i) 
r sin o e i ( uot +<t*') £( G 



r sin0 o e _i(wot+ * o) C(f ,t) 
,,t) 2- (1 + ro cos 6 )Z(C,t) 



(29) 
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FIG. 5: (Color online) The geometric phase as a function of R with different z. Here the coupling strength W = lOwo and 
wo = 1. The initial states of the open system are the same as in FigQ] 
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FIG. 6: (Color online) The information flow M as a function of R with different z. Here the coupling strength W = 10wo and 
Uo = 1. The initial states of the open system are the same as in FigQ] 

where 

f(C,r) = e^[cosh(-^) + i s inh(-^)], 

n = VI - ac. 

The parameters C and r are defined as C = and t = "ft. 

By the same procedure as in Sec lIII A} we establish relation between the geometric phase and the information flow, 
which is found to be similar to that in the last section. Hence, the conclusions for the relation between the geometric 
phase and the information flow hold true for open systems described by memory kernel master equation. See the 
numerical results shown in Figj8]and Fig|9l 

For the memory kernel master equation with exponential memory kernel, it is well known that, in some region of 
parameter C, the master equation Eq. (|26[) may lead to non-positive density matrix. In Ref . (23j . the positivity for a 
density matrix has been discussed: When C > j, the the memory kernel master equation with exponential memory 
is not valid, because the second perturbation used to drive the master equation does not suit this case. Recently, 
Breuer et al. has checked that this kind of master equation does not own memory effect when the positivity of the 
density matrix is conserved [24[ . When C < \ , there is no information flow back into the open system, no matter how 
to choose the parameter 70, so we choose 70 = O.Iluq and plot the geometric phase and the information flow Ai(C) in 
this region. Comparing Fig|8]with FigJH we can find the link between the information flow A4(C) and the geometric 
phase is the same as we found in Sec lIII Al 
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FIG. 7: (Color online) (a)The information flow M and (b) the geometric phase $gp as a function of R and the initial polar 
angle $o (in units of n). Here the coupling strength is W = IOcjo and wo = 1. The initial states of the open system are chosen 
as z — 1/2 and <po = n/6. 
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FIG. 8: (Color online) The geometric phase for the memory kernel master equation with exponential memory kernel as a 
function of the parameter C with different z. Here the dissipation constant takes 70 = O.Iljo and cjo = 1. We choose the same 
initial states as in Fig[T]for the open system to plot this figure. 
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FIG. 9: (Color online) The information flow M for memory kernel master equation with exponential memory kernel as a 
function of C with different z, where the dissipation constant is 70 = O.lwo and uuo = 1. The initial states of the open system 
are the same as in Fig[T] 
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IV. CONCLUSION AND DISCUSSIONS 



In summary, we have discussed the information flows and the geometric phase in different non-Markovian process. 
For an open two-level system with pure initial states, an analytic relation between the information flow and the 
geometric pahse has been given in terms of Bloch vector. For mixed initial states, two kinds of master equation, namely 
time-local master equation and memory kernel master equation with exponential memory, have been numerically 
studied. We find that, in both cases, the information flows influence the geometric phase directly, a relation between 
the geometric phase and information flow is numerically established. An understanding for the observation is provided. 

The forward and backward information flow are by definition different, but they complementarily describe the 
information exchange between the environment and the system. The backward information flow can be used to 
describe the Non-Markovianity of the open system, while the forward information flow was connected with the 
coherence loss. Neither forward information flow nor backward information flow can be measured directly, indicating 
that the measure of non-Markovianity defined in this way can not be directly observed in experiment. However, 
mixed state geometric phases are measurable and the feature caused by it has been observed [25j. This motivates 
the establishment of the connection between the geometric phase and the information flow. Indeed, The finding of 
this paper suggests that the geometric phase can reflect the non-Markovianity and then can serve as a measure of 
non-Markovianity for open systems. 

This work is supported by NSF of China under Grant Nos. 10775023, 10935010 and 10905007. 



APPENDIX 



In this appendix, we show in detail that the critical point is exactly the point where Af(t) begin to increase. By 
Eq.(fTTj). we write the integrand as 



A(t,R)=u [l+ - Tz ^ R) =], (30) 

2 2^/4D 2 {t,R)-2r z {t,R)-l l V ' 

where r z (t) = r(i) cos 6{t). To find the critical point, we take a derivative with respect to R 

j^A(t,R) = "^ 2 )3 go [(l + r cosg )| C (t,i?)| 2 + l]A| C ft,fi)|^ (31) 

noting that r(t,R) £ [0,1] and (1 + r a cos 9 a )\c(t, R)\ 2 + 1 > 1, we find, when -J^\c(t,R)\ 2 = 0, A(t,R) must reach its 
extremum. Moreover, we check the first derivative of the trace distance with respect to R, 

M D{t > R) = AD(t,R) [{1 + r ° cos0 o) 2 \ c (t, R)\ 2 + rl sin 2 6 ]-^\c(t, R)\ 2 . (32) 

Obviously, since D(t, R) £ [0, 1] and (1 + ro cos#o) 2 |c(i, R)\ 2 + r 2 , sin 2 9q > 0, we can see that the trace distance arrive 
at its extremum if and only if -^ \c(t, R)\ 2 =0, this is exactly the condition for the integrand to reach its maximum. 
Substituting Eq.© into Eq.([3"2|). we obtain 

^A^(i, i?) ~ ^A^(t, i?) = ^L-^[(l + r cose ) 2 |c(t, i?)| 2 + sin 2 ^o]^|c(*, ^)| 2 - (33) 

At the critical point, there is no backward information flow, i.e. Af(t) — 0. Thus if the critical point is the very point 
satisfied -J^\c(t, R)\ 2 =0, the derivative of A4(t) must be zero at this point. Because the property of the information 
flows, it is difficult to obtain an analytic results. The numerical result of Af(T) and M.{T) have been shown in FigflUl 
which validate our hypothesis. This is further confirmed by Fig. 1111 
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FIG. 10: (Color online) The measure A/" and M as a function of R with a pure initial state($o = 7r/3). The coupling strength 
W is chosen as 0.6. This figure shows that, at the start point (the point of crossover with the dot-dashed line) of the measure 
Af(T), the derivative of the measure M(T) with respect to R is zero. 
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